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A single cavity photon mode is expected to modify the Coulomb interaction of an electron system 
in the cavity. Here we investigate this phenomena in a parallel double quantum dot system. We 
explore properties of the closed system and the system after it has been opened up for electron 
transport. We show how results for both cases support the idea that the effective electron-electron 
interaction becomes more repulsive in the presence of a cavity photon field. This can be understood 
in terms of the cavity photons dressing the polarization terms in the effective mutual electron 
interaction leading to nontrivial delocalization or polarization of the charge in the double parallel 
dot potential. In addition, we find that the effective repulsion of the electrons can be reduced by 
quadrupolar collective oscillations excited by an external classical dipole electric field. 
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I. INTRODUCTION 

Experiments on atom or electronic systems of various 
types in photonic cavities have been opening up a new 
and exciting venue to test and manipulate photon-matter 
interactions in the strong coupling limit 

Modeling of these time-dependent systems has com¬ 
monly been aimed at evaluating their steady state prop¬ 
erties, that can include such diverse observables as their 
conduct ance,l^ the life tirne,!^ the broadening,!^ or the en¬ 
ergy shift of certain many-body states or modes of rel¬ 
evance. Here, we would like to draw attention to how 
the experimentally challenging transient time regime can 
reveal information about the interaction of the partic¬ 
ipating constituents, in our case, electrons and cavity 
photons. 

Many different formalisms have been used to describe 
transport in many-body systems on the nanometer scale. 
Most common have been m ethods built on various Green 
function approaches) ^^ * ^^ ! Les s common has been the 
use of master equations,!^^^ or a generalized master 
equation in combination with configuration interaction 
(Gi) [Ill where the Hamiltonian matr ix is diagonalized in 
a truncated Fock many-body space.!!^^ Between these 
two approaches is an interesting correspondence. When 
dealing with a system with two types of interactions, 
like the electron-electron interaction, and the electron- 
photon interaction, one typically may start by building 
a Green function for the electron including the mutual 
electron interaction. A second step would be to eval¬ 
uate the Green function of the full system dressed by 
the photons.!!^ Similarly, in a Cl-approach, one starts by 
diagonalizing the Hamiltonian matrix for the Coulomb 
interacting electrons, constructs a new many-body Fock- 
basis from the states of the Coulomb interacting electrons 
and possibly the photon number operator, and diagonal¬ 


izes then the full Hamiltonian matrix.!!^ This approach 
has been termed: “Stepwise introduction of model com¬ 
plexity”. Theoretically, similar structure and accuracy 
could be accomplished for the static properties in the 
two very different formalism, though their actual poten¬ 
tial depends strongly on the system of interest. Their 
difference becomes more obvious for the dynamic or the 
transport part. 

Here, our goal is to apply the latter strategy to ana¬ 
lyze the effective interaction of electrons in a very special 
system. We focus our attention on the interaction of two 
electrons in a parallel double dot system. We want to an¬ 
swer the question if it is possible to determine how the in¬ 
teraction of the two electrons with a single photon mode 
in the system modifies their repulsion. How is this effec¬ 
tive electron-electron interaction? Does it depend on how 
the question is approached, using static or dynamic prop¬ 
erties of the system. This question naturally arises when 
the energy spectrum of the Jaynes-Gummings modeP^is 
evaluated for strong electron-photon coupling and pho¬ 
ton energy out of resonance with the two atom levels.!^ 
In this case a strong negative dispersion is found with 
increasing atom-photon coupling, but the model was not 
intended for the off-resonance situation, and it is known 
that the diamagnetic electron-photon interaction has to 
be included in a related model of ele ctrons and cavity 
photons in a short quantum wire.!^^^^ 

We explore these question for the case of a parallel dou¬ 
ble dot system with two interacting electrons closed with 
respect to electron exchange with the environment, and 
compare the results with the situation in which the sys¬ 
tem is the open version of the model, when it is coupled 
to two external leads with a bias voltage that promotes 
current through it. We hope to convince the reader that 
the geometry of this particular system furnishes it with 
very special properties worth studying. 
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The paper is organized as follows: First we introduce 
the components of the static model, the many-body for¬ 
malism for exactly interacting electrons and cavity pho¬ 
tons in a confined system with a continuous geometry 
describing two parallel dots embedded in a short quan¬ 
tum wire. We explore its static and dynamics proper¬ 
ties with increasing electron-photon interaction. The dy¬ 
namic properties are obtained by observing the response 
of the system to an excitation with an external classi¬ 
cal electric dipole pulse of short duration. We then pro¬ 
ceed to open the system up to transport by coupling it 
to two external electron leads and investigate the effects 
of varying the electron-photon interaction on the time- 
dependent currents to and from the system, the charge 
in it, and its distribution through it. 


II. CLOSED SYSTEM 



In the closed version of the system we consider two 
electrons in the parallel double-dot potential landscape 
shown in Fig. interacting with the Coulomb interac¬ 
tion that is accounted for via a configuration interaction 
(Cl) approach in the basis of noninteracting two-electron 
states.l^ We assume GaAs parameters, k, = 12.4, and 
m* = O.OGTmg. 

The potential landscape of the short quantum wire and 
the embedded dots is described by 


V(x,y) 
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with MIq = 2.0 meV, Vd = —6.5 meV, [3 = 0.03 nm“^, 
di = —50 nm, ^2 = —50 nm. Lx = 150 nm, and 0 is the 
Heaviside step function. 

The Hamilton operator for the closed system is 




= j d?rxl)^ (r) I + ^(r) I V’(r) + Hem + -f^Coui 
-IJ d0 j(r) ■ j d\ p{v)A0 (2) 


FIG. 1. Schema of the leads, the central system with double 
dots and the cavity (upper panel). The transport direction 
is indicated by magenta arrows, and the photon mode with 
golden arrows. The potential landscape defining the two par¬ 
allel quantum dots in a short parabolically confined quantum 
wire (lower panel). The effective magnetic length — 23.8 
nm, hLLo — 2.0 meV, and H = 0.1 T. 


our results. i^EM is the Hamiltonian for the cavity pho¬ 
ton mode with energy and the electron-electron in¬ 
teraction represented by i^coui- The vector potential of 
the cavity photon field is A^. The two last terms of 
the Hamiltonian © are the para- and the diamagnetic 
electron-photon interactions, respectively. The charge 
and the charge-current densities are 

p=-eip^'ip, j =, (4) 

The natural length scale defined by the parabolic con- 
fi nement an d the external magnetic field is = 
\/h /with the cyclotron frequency ujn = 
eH/(m*c), and ftl = uj^ ^ ftl 

We assume a rectangular photon cavity (x^y^z) G 
{[—ac/2, ac/2] x [—ac/2,ac/2] x [—dc/2, (ic/2]} with the 
short quantum wire centered in the z = 0 plane. In 
the Goulomb gauge the polarization of the electric field 
is chosen parallel to the transport in the x-direction by 
selecting the TEqh mode, or perpendicular to it by se¬ 
lecting the TEioi mode 


with 


TT = (p + ^Aext) , (3) 

where Aext is an external vector potential producing a 
homogeneous small magnetic field H = 0.1 T along the 
2 :-axis included in the calculations in order to break spin 
and orbital degeneracies to guarantee accuracy of the re¬ 
sults. We have tested that this magnetic field can be 
made vanishingly small without qualitatively changing 


A^(r) 



.4- -j- 



cos 



( 5 ) 


for the TEoii and TEioi mode, respectfully. The 
strength of the vector potential, A, determines the cou¬ 
pling constant ^em = eALlwdw/c leaving a dimensionless 
coupling or polarization tensor 


gij = ^{{i\e-7T\j)+h.c.}, (6) 
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where \i) is a single-electron state of the short quantum 
wire. We are describing an electron system with many 
electron states that will generally not be in resonance 
with the cavity photon mode so we retain all resonant 
and antiresonant terms in the photon creation and an¬ 
nihilation operators and a (not applying the rotating 
wave approximation). We assume the cavity to be much 
larger than the electron system and neglect the spatial 
variation of the vector potential A^. In a Fock basis for 
the Coulomb interacting electrons we arrive at the 

Hamiltonian 

Hs = '^ + tiLoa^a 

+ gEM gij {a + a^} 

^vij 

^ ^ ^ livi 

I - (aa + a^a^) | , (7) 

where di and dj are the electron annihilation and creation 
operators for the single-electron state |i), and the Fock 
basis for noninteracting electrons is connected to 

the interacting electron Fock basis by the unitary trans¬ 
form l/i) = The spectrum and states of this 

Hamiltonian Q are found by diagonalization in a Fock 
space constructed by a tensor product of the eigenstates 
of the photon number operator a and the Fock 

basis of Coulomb interacting electrons In this basis 

without electron-photon interactions each electron state 
can be found with a different number of photons, these 
state are called photon replicas of the electron states. 
The electron-photon interaction makes the photon con¬ 
tent of an interacting state deviate from integer numbers 
for strong coupling. If the deviation is not large the “pho¬ 
ton replica” concept can still be used for the interacting 
system. The “first photon replica” of a certain state is 
then the tensor product of that state and the eigenstate 
of the photon number operator with one photon. The di- 
agonalizations for the Coulomb interaction on one hand, 
and for the electron-photon paramagnetic- and diamag¬ 
netic interactions are done separate ly in truncated bases 
to optimize the size needed. 

Here, we come to an essential point in our approach: 
Researchers have found that for the Jaynes-Camming 
modeP^ in the strong coupling limit it can be difficult to 
achieve convergence in a basis constructed with the pho¬ 
ton number operator.!^ We do not encounter this prob¬ 
lem for our model when the diamagnetic electron-photon 
interaction is included ([^,^ but as the cavity-photon 
dressed electron-electron interaction polarizes the charge 
density we have to be very careful about including enough 
of higher energy electron states in our bases. 

The polarization of the photon field will be chosen 
along the x-axis parallel to the electron transport when 
the system will be opened, or perpendicular to it, along 


the ^-axis. When calculating matrix elements we assume 
the wavelength of the cavity photons to be much larger 
than the size of the electronic system. At the moment 
the system is closed with respect to electrons and pho¬ 
tons. Since, later we intend to open up the system we 
include in our calculation states with up to 3 electrons. 
The energy spectra for the 60 lowest states together with 
their photon and spin content are displayed in Fig. [^for 
the two polarizations of the photon field. 



10 20 30 40 50 60 
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FIG. 2. The energy spectra for the 60 lowest many-body 
states /i = 1, 2, • • • , 60) showing their electron (e), pho¬ 
ton (7), and spin content (s^). The photon field is x-polarized 
in the upper panel, but y-polarized in the lower. The left and 
right chemical potentials, fiL and /j.r, are shown as horizon¬ 
tal lines for the case of the open system. HQq = 2.0 meV, 
hu = 0.8 meV, ^em = 0.1 meV, B — 0.1 T, = 12.4, and 
m* = 0.067me. 

The inclusion of the electron-photon interaction causes 
dispersion of the energy for both one- and two-electron 
states with respect to the electron-photon coupling con¬ 
stant ^EM-^The former can be thought of in terms of an 
effective potential, and the latter we explain as an effec¬ 
tive cavity-photon dressed electron-electron interaction. 

We start our study by analyzing the spatial extent of 
the charge density of the many-body relevant states. In 
the upper panel of Fig.j^we present the average extent of 
the first photon replica of the 2-electron ground state as 
a function of the coupling ^em for the closed system with 
the photon energy huu = 0.8 meV, compared to the case of 
hw = 2.0 meV when the ^-polarized photon field is in res¬ 
onance with the confinement energy of the system in the 
^-direction. The dressing of the electron-electron interac¬ 
tion by the paramagnetic electron-photon term stretches 
the charge density in the direction of the polarization of 
the photon field while the dressing by the diamagnetic 
term spreads it in all directions. Here, we see that in all 
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FIG. 3. (Upper panel) The expectation value {{r/aw)‘^) for 
the closed system at t = 0 for x- and y-polarization of the 
photon field. (Lower panel) the expectation values {{x/aw)‘^) 
and ((y/a^(;)^) for the closed system at t = 0 for x- and y- 
polarization of the photon field. The photon energy is Uem = 
hw = O.S meV except for one curve for comparison in the 
upper panel. HQq = 2.0 meV, and B — ^.1 T. 


cases the extent of the charge density increases, though 
by a small amount for the photon field out of resonance. 
In the lower panel of Fig. [^we analyze the change in the 
extent of the charge density of the same state, but along 
the X- and ^-axes. Surprisingly, the charge density of the 
2-electron state always shrinks in the x-direction and ex¬ 
pands in the ^-direction with increasing ^em, irrespective 
of the polarization of the photon field. 

To show how this happens for the x-polarized photon 
field we compare in Fig. I^the probability density for the 
first photon replica of the 2-electron ground state with 
the changes in the density defined as 5n = n(^EM = 
0.16) — n(0). Clearly, the central section of the density 
shrinks in x-extent, but together with the stretching in 
the ^-direction we notice a stretching in the x-direction 
for the density located in the quantum dots. 

From these observations we conclude that the polariza¬ 
tion of the charge or probability density of the 2-electron 
state indicates an enhanced repulsion between the elec¬ 
trons. 

Yet, another way to gauge the effective Coulomb inter¬ 



FIG. 4. (Left panel) The probability density for ground state 
for 2 electrons and 1 photon with ^em = 0. (Right panel) The 
difference in the charge density 5n — n(^EM = 0.16) — n(0) for 
the same state as in the left panel. a:-polarized photon field, 
huj = 0.8 meV, B = 0.1 T. 


action is to excite the system by a weak classical electrical 
dipole field polarized in either direction 

Kxt(r, t) =Vt{r- e) exp {-Tt) sin {ojit), (8) 

with Vt = 0.1 meV, hV = 0.1 meV, and huji = 0.4 meV, 
and e determines the polarization of the pulse. We are 
interested in the time evolution of the system after it 
is subjected to a short dipole excitation pulse that can 
be polarized in the x- or ^-direction and adds only a 
tiny amount of energy to it.^^ Since our system has a 
nontrivial geometry and interactions between the elec¬ 
trons compared to the circular dot explored in an ear¬ 
lier publication with only Coulomb interactiorP^we know 
that the system will respond not only with dipole oscilla¬ 
tions. There are higher order modes and Fig. shows the 
Fourier transform of the expectation values for x‘^ and y‘^ 
for ^-polarization of the photon and the classical exci¬ 
tation field. In both cases the fundamental modes (de¬ 
picted in Fig. show a weak softening with increasing 
coupling ^^em- Since we are using a Cl-formalism for the 
interaction we can identify which many-body states par¬ 
ticipate in the oscillations, their occupation being shown 
in the left panel of Fig. For low ^fEM it is mainly |61), 
a state with peaks in the probability density located in 
the dots, but pushed a bit to higher values of y than in 
116) and then a considerably more density in the center 
between the dots, see Fig. Comparing the probability 
densities in the left panel of Fig. and we see that the 
weak dipole excitation of the system polarized in the y- 
direction causes internal motion of the charge or density 
with the charge in each dot oscillating in antiphase. 

Furthermore, and again due to the Cl-formalism, when 
we compare the dispersion for the energy of the states 
|16) and |61) in the right panel of Fig. we see that 
the energy difference of the two states decreases with 
increasing ^em reflecting the softening of the modes in 
Fig.m From this analysis of the closed system excited by 
a weak classical dipole electric field we conclude that the 
effective Coulomb interaction is reduced by the photon 
field, contrary to our conclusion earlier when no energy 
was added to the system. 
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FIG. 6. (Left) In the closed system (after the excitation pulse 
vanishes) the occupation of the states as a function of the 
electron-photon coupling ^em- (Right) The dispersion of the 
initial state |16) and the excited state |61) (both with 2 elec¬ 
trons and 1 photon), y-polarization of the photon field and 
dipole excitation pulse in the same direction. States |8) and 
129) are 2-electron states with 0 and 2 photons, respectively. 
Eem = huj = O.S meV, hQo = 2.0 meV, and B = 0.1 T. The 
arrows indicate the respective equidistant y-axis. 


FIG. 5. The Fourier transform of the expectation values 
{(x/awy){t) (upper), and (lower) as a function 

of the electron-photon coupling ^em • y-polarized photon field 
and excitation pulse. hQo = 2.0 meV, and B = 0.1 T. 

III. OPEN SYSTEM 

We now turn our attention to the open system. To 
open up the system to semi-infinite left and right leads 
with parabolic confinement in the ^-direction we intro¬ 
duce a time-dependent coupling iLy( t) bet ween the cen¬ 
tral double-dot system and the leads 

HTit) = + {Tl,rdlc,i} , (9) 

i,l ^ 

where di annihilates an electron in the single-electron 
state of the the central system labelled with i, and 
creates one in lead I with momentum and subband in¬ 
dex labelled by the composite index The coupling 
Tg- depends on the shape of the wavefunctions for each 
single-electron state in the leads and the central system in 
the contact areas defined between the subsystems. The 
time-dependence of the coupling is determined by the 
switching functioiP^ 



for both leads with a = 0.3 ps“^. The state dependence 
of the coupling tensor in the single-electron picture is 
carried on to the corresponding coupling tensors in the 
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FIG. 7. The probability density for 2 electrons and 1 photon 
with ^em = 0.05 meV in the excited state |61). y-polarized 
photon field, HQq = 2.0 meV, and B = 0.1 T. 

many-body Fock spaces used in the calculation.E^l^ The 
Liouville-von Neuman equation for the density operator 
of the total system is then projected on the centra l syste m 
by tracing out variables pertaining to the leadsIn 
this proceedure we extract an overall coupling constant 
= 0.5 meV. In the generalized master equation 
(GME) for the reduced density operator 

ps(t) =-^[Hs,Ps{t)]-l dt'K[t,t-t';ps{t')] (11) 

we keep up to second order terms of the lead-system cou¬ 
pling Ht in the integral kernel of the dissipation term 

K[t,s-,ps{t')]=Tr^n{[HT{t), [U{s)HRt')U+(s), 

C/s(s)ps(i')G(s)/5L/9R]]}, (12) 
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where Us{t) = exjp {—iHst/h} is the unitary time evolu¬ 
tion operator for the closed system, and U (t) is the cor¬ 
responding operator for the uncoupled system and the 
leads. TrLR is the trace operator for variables of the 
left and right leads, and pl and pl are the density ma¬ 
trices for the leads before coupling to the central sys¬ 
tem. The time-dependent reduced density operator can 
be used to ca lculat e mean values of measurables of the 
central systenf^^^ under the influence of the leads, that 
serve as electron reservoirs. Often, only properties of 
the steady state are of interest to evaluate the conduc¬ 
tance of the system, life-time of certain states, their shift 
or broadening. Here, we are interested in analyzing the 
interaction of electrons that is best explored in the tran¬ 
sient regime where the non-Markovian character of the 
GME is still significant. 

In the numerical calculations for the closed system we 
commonly need the size of various bases to be on the or¬ 
der of few thousand states. Due to the complex structure 
of the dissipation terms in the GME we have to limit the 
size of the final many-body basis used for the transport 
calculation to 120-192 states. The time-evolution of the 
closed system up to t = 1200 ps for one set of parameters 
needs 1-2 days on 16 GPU and 512 GPU threads, and for 
the open system up to 300 ps it needs 5-7 days on 12- 
16 GPU threads. The parallelization of the dissipation 
terms of the GME is a difficult task and could definitely 
be improved. All programs are written in EORTRAN 
2008 with OpenMP to maximize their speed. 

We begin the exploration of the open system by select¬ 
ing the bias window to include the first photon repiica 
(116)) of the two-eiectron ground state (|8)). The pho¬ 
ton energy is hw = O.S meV, not ciose to a resonance for 
any iow lying 2-electron states. Inspection of Eig.j^shows 
that this photon energy is not close to a resonance for any 
one-electron state for the case of x-polarization, but close 
to a resonance for the one-electron ground state |1) bound 
in the dots and the one-electron state |3) just above the 
dots in the case of ^-polarization. This can be seen by 
the fact that the mean photon number for |3) and |5) and 
their opposite spin counterparts is not close to an inte¬ 
ger, reflecting a Rabi-splitting and entanglement.!^ The 
resulting extra polarizability of the one-electron states 
for the ^-polarized photon field is not seen for the low 
energy 2-electron states. Eigure shows the mean elec¬ 
tron number in the system as a function of time. In 
Eig. I^c) and (d) we observe how one of the electrons 
leaves the system when initially the first photon replica 
of the two-electron ground state is occupied. In panels 
(a) and (b) we see how a second electron enters the sys¬ 
tem when initially the first photon replica (|6)) of the 
one-electron ground state (|1)) is occupied. The results 
in Eig. show that for the x-polarized photon field the 
depopulation of the 2-electron states is enhanced with 
increasing electron-photon coupling and the addi¬ 

tion of a second electron is resisted by increased ^em- 
The trend is not linear with ^em- Similar trend is seen 
for the ^-polarization, but there the nonlinearity is in- 



0 100 200 300 0 100 200 300 
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FIG. 8. The average electron number (Ae(t)) as a function of 
time in an open system in the case of one electron initially in 
the first photon replica of the one-electron ground state in the 
system for x-polarized photon field (a), y-polarized photon 
field (b), and for two electrons initially in the first photon 
replica of the two-electron ground state in the system for the 
case of x-polarized photon field (c), y-polarized photon field 
(d). iiL = 1.4 meV, /xr = 1.1 meV, ghRaH^ = 0.5 meV, 
hQo = 2.0 meV, Eem — fvuo — 0.8 meV, and R = 0.1 T. 


creased, and for some regime of low coupling it appears 
to be easier to add the second electron into the system. 

The conclusion is that in this situation (marked by the 
location of the bias with respect to the energy spectrum 
and the geometry of the system) the repulsive Goulomb 
interaction is enhanced by the dressing cavity photon 
field. An analysis of the polarization of the charge den¬ 
sities by the photon field gives supporting arguments 
for this behavior. The dressing by the cavity-photon 
field changes the shape of the wavefunctions and the 
charge density in the contact areas to the leads ,1^ but 
due to the intricate interplay of the dressing by the para- 
and dian iagnetic electron-photon terms with the system 
geometr}E^l^^l^ the polarization of the charge densities 
is nontrivial, and in addition one has to account for the 
changing occupation of the one- and two-electron states. 

We have verified that qualitatively similar results are 
obtained for the two photon replica of the 2-electron 
ground state. The same can be said about the 2-electron 
ground state except there {x^) decreases slightly with in¬ 
creasing ^em, and {y‘^) only slightly increases. Eor the 2- 
electron ground state the dispersion of the Eourier trans¬ 
form of {x^) and {y‘^) is fairly constant (corresponding 
to Eig. 1^, while the transport (corresponding to Eig. 
indicates enhanced repulsion. 


IV. SUMMARY AND DISCUSSION 

Summarizing, we emphasize that this analysis is only 
possible by taking account of both the para- and diamag- 











7 


netic electron-photon interactions and the geometry of 
the system through the wavefunctions used to construct 
the various bases used in the calculationThe parallel 
double dot system brings a fine-tuned interplay between 
the para- and diamagnetic interactions that is not seen 
to the same extent in systems with a single dot or se¬ 
rial dotsp^ even though both interactions are important 
to describe phenomena when the photon field is not in 
resonance with the electron systemP^ To make this abso¬ 
lutely clear: Our results are model dependent and point 
out a particular system with interesting properties. 

We notice that the answer to the question about the 
infiuence of the cavity photon mode on the effective 
electron-electron interaction depends on the stand point 
we take to gauge it, or define it. The entity, two elec¬ 
trons dressed by cavity photons, appears as a quasipar¬ 
ticle with different properties in different “setups”. 

In our parallel double dot system investigated here the 
dispersion of the energy of two-electron states bound to 
the dots shows that the repulsion is enhanced by the 
electron-photon interaction. This is not a universal be¬ 
havior as can be seen for two electrons in a quantum wire 
without the embedded dots.l^ We conclude thus that the 
effects of the photon field in the cavity could not gener¬ 
ally be described by a constant renormalization of the 
dielectric constant of the surrounding material. 

We find that in our transport setup the photon field 
enhances the Coulomb repulsion of the two electrons in a 
parallel double quantum dot. This is true for initial states 
below, in, or above the bias window, as the coupling to 
the leads at t = 0 and the photon field enables slight cou¬ 
pling to states higher in the energy spectrum, and partic¬ 
ularly to a one-electron state lower in the spectrum, i.e. 
the one-photon replica of the one-electron ground state. 

In the closed system, in equilibrium, when gauged 
through the extent of the charge densities of the two- 
electron states, the effective interaction of the two elec¬ 
trons is increased, or left untouched. On the other hand, 
when probed by observing the dispersion of the internal 
collective oscillation modes, excited by a weak classical 
electric dipole pulse, the effective interaction of the two 
electrons is reduced, or left unaffected. 

When we observe the strengthening of the repulsion 
between the two electrons due to the dressing of their 
mutual interaction by the cavity photon field it is seen 
by the spreading of the charge (or polarization), but in 
the case of the transient transport we have to bear in 
mind the infiuence of occupation of states with one more 
or less electron than the initial state had. 

A clear exception to the enhancement of the repulsive 
interaction by the electron photon coupling in the present 
model is the softening of the relative component of the 
collective oscillation mode as measured by the Fourier 
transform of the expectation values of and in Fig. 

Here, Fig. compared to Fig. indeed shows that 
in addition to a spreading of the charge density both in 
X- and ^-direction we see a very small accumulation of 
charge between the quantum dots during the oscillations. 


This is valid for both linear polarization of the photon 
field in the x- and ^-direction. The charge densities show 
a quadrupole oscillation that is facilitated by the two 
terms of the electron-photon coupling thus resulting in a 
softening of the mode with higher coupling. 

The Fourier transform of the dipole response as mea¬ 
sured by {x) and (y) is much less affected by the coupling 
giving weight to our argument that the softening of the 
relative mode can be referred to the interaction between 
the electrons. 

One can wonder if the changes in the effective electron- 
electron interaction observed here can be attributed to 
the Pauli principle and would thus also be seen in non¬ 
interacting systems. We are sure the electron-photon in¬ 
teraction would have effects on electrons that were not 
Coulomb interacting, but that could be explained as the 
emergence of an electron-electron interaction. Generally, 
it is impossible to distinguish effects caused by changes 
in the structure of the state space and the underlying 
wavefunctions. One might even ask if such an attempt 
would depend on the basis used. An increasing complex¬ 
ity of a state space can be used to describe the geometry 
of a system. 

The system we explore here is a complex system with 
many parameters requiring long CPU-time on paral¬ 
lel facilities. We find dispersion of the energy of the 
two-electron states indicating increased repulsion with 
stronger electron-photon coupling. As mentioned before 
this is not always the case for two-electron systems in a 
photon cavity.^ We are not only gauging these systems 
in their equilibrium and it is likely that variants with 
different geometry can sustain different collective modes 
leading to a reduction of the effective electron-electron 
interaction as we have found here. 

It is interesting to notice that at two different points 
in the model we are witnessing an effective light-light 
scattering or interactions. First, through the dressing 
of the electrons by the cavity photons as the Coulomb 
interaction of the electrons is also mediated by photons. 
Second, the diamagnetic electron-photon interaction nec¬ 
essary here is essentially a nonrelativistic light-light in¬ 
teraction that in solid state systems can lead to nonlinear 
effects for much lower strength than is needed in a vac¬ 
uum treated relativistically. 

The steady state in the open system only depends on 
the location of the bias window in the energy spectrum 
of the system and the coupling strength. It would thus 
be independent of the initial number of photons in the 
cavity. Our results for the transient regime depend on 
the initial number of photons and electrons in addition to 
the location of the bias window, though we have verified 
that our finding here is valid for an initial state with few 
photons (0-2). We are describing a system with electron 
states, coupled to photons, that can have very different 
coupling to the leads due to their geometry making the 
estimation of the time needed to reach a steady state 
difficult. This together with the fact that conceivably 
some quantum electro-optical devices might be operated 


in the late transient time in order to save time and reduce 
dephasing lets us focus on that regime. 
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